Testing "microscopic" theories of glass-forming liquids 
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We assess the validity of "microscopic" approaches of glass-forming liquids based on the sole 
knowledge of the static pair density correlations. To do so we apply them to a benchmark pro- 
vided by two liquid models that share very similar static pair density correlation functions while 
displaying distinct temperature evolutions of their relaxation times. We find that the approaches 
are unsuccessful in describing the difference in the dynamical behavior of the two models. Our 
study is not exhaustive, and we have not tested the effect of adding corrections by including for in- 
stance three-body density correlations. Yet, our results appear strong enough to challenge the claim 
that the slowdown of relaxation in glass-forming liquids, for which it is well established that the 
changes of the static structure factor with temperature are small, can be explained by "microscopic" 
approaches only requiring the static pair density correlations as nontrivial input. 
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I. INTRODUCTION 

Despite decades of intensive research, the connection 
between dynamics and structure in glass-forming liquids 
remains elusive. As a matter of fact, this is likely one 
main reason that explains the difficulty in producing a 
comprehensive theory of the glass transition 1 -. The struc- 
ture of glass-forming liquids, as experimentally probed 
through the static structure factors or equivalently the 
pair correlation functions, does not appear to change 
much when temperature is decreased toward the glass 
transition while under the same condition, the relax- 
ation time dramatically increases by orders of magni- 
tudes. To incorporate this, rather central, observation, 
several working hypotheses have been proposed: (i) do 
away with the statics and assign the slowdown of re- 
laxation to purely dynamical features such as emerg- 
ing kinetic constraints^" 3 -, (ii) focus on structural quan- 
tities that go beyond the pair density correlations and 
represent higher-order and subtler correlations between 
particles 4 - - — , (iii) put forward mechanisms by which the 
small modifications of the static pair density correlations 
can be extremely amplified to produce spectacular dy- 
namical changes. The theories based on (i) and (ii), 
such as those involving dynamical facilitation 3 -, structural 
frustration—*^ or a random first-order transition 4 * 1 ^, 
usually require some sort of phenomenological input, ei- 
ther because they involve some coarse-graining that in- 
troduces effective parameters only indirectly related to 
microscopic details or because they are intractactable in 
their full-blown generality and need approximate treat- 
ments. The paradigm of a type (iii) theory is the mode- 
coupling theory of the glass transition 14 that predicts a 
freezing of the dynamics through a nonlinear feedback 
effect affecting the fluctuations of density. 

The appeal of so-called "microscopic" theories is that 
they are able to make predictions from the knowledge 



of the molecular properties, among which the interac- 
tion potentials 1 ^. In the case of glass formation where 
it seems likely that some form of collective or cooper- 
ative behavior is at work and, correlatively, that some 
form of independence from the molecular details charac- 
terizes the phenomenology, it is not clear that this should 
necessarily be the best route to follow. In addition, "mi- 
croscopic" theories of the glass transition involve strong, 
uncontrolled approximations that are usually not justi- 
fied at a microscopic level. Nonetheless, such theories 
are at least amenable to crisp tests on atomistic model 
systems. 

In practice, "microscopic" approaches of the dynamics 
of glass-forming liquids are built on the knowledge of the 
static two-body density correlation functions (the static 
structure factors in Fourier space). The latter encode 
the microscopic information but represent only a partial 
description of the liquid structure. The issue we address 
here then boils down to the following question: how far 
can one go with the idea that the small observed changes 
in the static structure factors are sufficient to describe 
the slowdown of relaxation in glass-forming liquids? 

The approaches we consider in this paper, with var- 
ious degrees of exhaustivity and depth, are the mode- 
coupling theory (MCT), for which we shall mostly sum- 
marize the results obtained in our previous study 16 , 
the "microscopic" implementations of the random first- 
order transition (RFOT) theory through density func- 
tional theor y 17 ' 18 or replica formalis m 19 ' 20 , an approach 
that mix aspects of the MCT with thermally acti- 
vated events^ 1 -, and finally a proposed correlation be- 
tween the dynamics of a glassformer and its static 
pair correlation functions via the two-body contribu- 
tion to the excess entropy 22 . To test these approaches, 
we have used as a benchmark the glass-forming sys- 
tems that we have already extensively investigated by 
computer simulatio n 23 ' 24 : the Kob-Andersen binary 
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Lennard- Jones model and its WCA (for Weeks- Chandler- 
Andersen 25 ) reduction to truncated, purely repulsive po- 
tentials. These two models have been shown to share very 
similar static pair density correlation functions while dis- 
playing increasingly distinct evolution of their relaxation 
times as temperature decreases. A wide range of densi- 
ties can moreover be considered, from low densities near 
the liquid-gas spinodal of the full binary Lennard- Jones 
mixture to high densities at which the difference between 
the two models progressively vanishes, through the con- 
ventional liquid/supercooled liquid range. 

We find that the approaches based on the pair corre- 
lation functions, or equivalent ly on the static structure 
factors, are unsuccessful in describing the difference in 
the dynamical behavior of the two glass-forming models, 
the binary Lennard- Jones mixture and its WCA reduc- 
tion. This is true even after applying for some theories a 
global rescaling of the predictions to better fit the simu- 
lation data for one model, say the binary Lennard- Jones 
one (in absolute terms, the predictions are indeed never 
very good) : the difference in the slowdown of relaxation 
of the two models is still not correctly accounted for. As 
already mentioned, our study is not exhaustive, and we 
have not tested the effect of adding corrections to the de- 
scription in terms of static pair correlations by including 
for instance three-body density correlations. In addition, 
we have chosen to focus on simple glass-formers. Yet, 
our results appear strong enough to challenge the claim 
that the viscous slowing down of glass-forming liquids, for 
which it is well established that the changes of the static 
structure factor with temperature are small, can be ex- 
plained by "microscopic" approaches essentially based on 
the knowledge of the static pair density correlations. 



II. BENCHMARK SYSTEM 

The glass-forming systems that we consider are the 
three-dimensional Kob-Andersen binary Lennard- Jones 
mixture 26 (denoted LJ in the following) and its reduc- 
tion to the purely repulsive part of the pair potentials 
proposed by Weeks, Chandler and Andersen^ 5 - (denoted 
WCA in the following). These are 80:20 mixtures of A:B 
atoms with interatomic pair potentials 
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i a g is equal to the position of the 
minimum of v a fj(r) for the WCA model and to a con- 
ventional cutoff of 2.ha a p (merely introduced for prac- 
tical reasons with no impact on the physical quantities) 
for the standard LJ model; C Q( g is a constant that is 
fixed such that ^(r^) = 0. The Molecular Dynamics 
simulations have been performed in the NVE ensem- 
ble, after equilibration at the chosen temperature, with 
N = 1000 particles, and we have studied a broad range 



of densities p from 1.1 to 1.8 (a detailed description of 
the phase diagram is given in Refi^). Lengths, tempera- 
tures, and times are given in units of a a a, ^AA/ksi and 
{mo\ A /4&tAA) 1 l 2 , respectively. 

These systems form a benchmark for the kind of inves- 
tigation that we want to carry out as their equilibrium 
pair structure is very close while their dynamics strongly 
diverge as temperature is lowere d 23 ' 24 . Note that the 
two models should be considered at the same density 
(pressure is very sensitive to the presence or absence 
of attractive forces). The observables that we measure 
in the simulation are the partial static structure factors 
S a p(q) and pair correlation functions g a f}{r) as well as 
the self-intermediate scattering functions F"(q,t) with q 
near the position of the peak of the static structure fac- 
tor at the most commonly studied liquid density p = 1.2 
(q<JAA — 7.2); we extract the relaxation time from the 
latter, with the conventional choice F"(q,t = r) = 1/e. 



III. TESTED "MICROSCOPIC" APPROACHES 

The "microscopic" approaches of the dynamics of 
glass-forming liquids that we assess by comparing to sim- 
ulation data on the above benchmark systems comprise 
first what can be taken as bona fide theories, the mode- 
coupling theory (MCT)^ 4 - and the random first-order 
transition (RFOT) theoryi 3 .. 

The MCT (for more detail, see our previous article^ 6 -) 
is based on a nonlinear differential equation for the time 
evolution of the two-point correlator of the density fluc- 
tuations. It can be derived for a Newtonian or for a 
Brownian dynamics, and, since the latter is somewhat 
simpler and leads to the same behavior at long times, 
this is the one that we have used. The time dependence 
of the intermediate scattering functions (two-point cor- 
relator of the density fluctuations) for the liquid mixture 
is then governed by 



D q 2 S- 1 (q)F(q,t) 



-J dt'M{q,t-t')^F{q,t'), 
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where Dq is the diffusion coefficient of an isolated Brow- 
nian particle and F(q, t) is the matrix formed by the 
intermediate scattering functions F a p(q,t) [such that 
Fa/3{q,t — 0) = S a p(q)]. The physics is contained in 
the so-called memory kernel M(g, t) whose explicit form 
is derived through a series of approximations from the 
exact but intractable formal expression; it is a functional 
of the intermediate scattering functions F(q,t) (taken at 
the same time t) and, in the common implementation of 
the theory, of the static structure factors S(q): 
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where the expression of M^^^ 1 5 i s given in Ref. [27| . 
(Note that in principle this expression involves the static 
triplet correlation function but the latter is almost al- 
ways considered in a factorized approximation.) In this 
formulation, the only nontrivial input is therefore the 
partial static structure factors S a p{q)- An equation can 
also be derived for the intermediate scattering functions 
F"(q,t)2&. However, it is known that for wave- vectors 
corresponding to typical interatomic distances the tem- 
perature dependences of the relaxation times associated 
with the collective and the self-intermediate scattering 
functions are similar and, moreover, that the predicted 
critical temperatures T c at which the dynamics freezes 
are identical. 

The second type of theory that we shall test are micro- 
scopic implementations of the random first-order transi- 
tion (RFOT) theory. The latter builds on an analogy 
with the behavior of mean-field spin glasses without re- 
flection symmetr y 29 i 30 . It provides a scenario for the 
slowdown of relaxation in glass-forming liquids that fo- 
cuses on the free-energy landscape, which is postulated 
to be characterized by an exponentially large number of 
"metastable" states and the associated configurational 
entropy (or complexity), and on an entropy-driven ac- 
tivated relaxation, which leads to a heterogeneous "mo- 
saic" liquid stated. Microscopic implementations of the 
RFOT theory can be realized at a mean-field level where 
free-energy "metastable" states have an infinite life-time 
and are therefore properly defined. They can be formu- 
lated either in terms of a density functional approach 
(DFT)iL±&2i. or w ithin a replica formalisinl 9 .! 2 ^ 3 ^ 3 . In 
these mean-field-like approaches, one does not have a di- 
rect access to the dynamics of the system but indirect 
information is provided through the configurational en- 
tropy and the two critical temperatures, Td at which the 
complexity jumps to a finite value and Tk at which the 
complexity vanishes. 

In the DFT treatment, one looks for the density pro- 
files that minimize some appropriately derived (mean- 
field) density functional. The main steps are the formu- 
lation of the density functional itself and the a priori 
characterization of the structure of the metastable states 
in terms of trial (amorphous or "aperiodic" 34 ) density 
profiles. The commonly considered Helmoltz free-energy 
functional is the Ramakrishanan-Yussouff 35 approxima- 
tion which involves a Taylor series expansion (truncated 
after the quadratic term) of the excess free energy around 
the liquid phase with mean density pq. For a one- 
component liquid, it is expressed as 

F[p] = J dvp{v) (ln[A 3 p(r)] - l) + F ex {p Q ) 
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where the first term of the right-hand side is the ideal-gas 
functional, with A the thermal wavelength, and J-" e:c (po) 
is the excess free energy of a homogeneous liquid of den- 



sity po; c(r; po) is the direct correlation function of the 
homogeneous liquid and is related to the pair correlation 
function g(r; po) through the Ornstcin-Zernike integral 
equation 36 . The density profiles p(r) of the metastable 
states are generally described through a sum of Gaussian 
functions centered about a given amorphous (or aperi- 
odic) lattice, with the width of the Gaussians character- 
ized by a variational localization parameter on which the 
minimization of the free energy is then performed. In 
practice, the amorphous lattice is chosen from a random 
close packing of hard spheres with some effective diam- 
eter. Additional approximations can also be introduced 
in order to further simplify the computation a 1 7 ' 18 ' 31 . 

An alternative approach uses the replica formalism, 
which allows one to directly study the statistical prop- 
erties of the metastable states and the configurational 
entropy. To this end, one introduces m replicas of the 
original liquid system and, if one is only interested in the 
liquid above the putative RFOT to an ideal-glass phase, 
one considers the difference between the free energy in 
the limit m — > 1 + (after having taken the thermody- 
namic limit) and the homogeneous liquid value obtained 
for m = l 19 ' 32 . At low enough temperature, one ex- 
pects the atoms in the m copies to cluster around around 
common positions in space, thereby forming a "molecu- 
lar bound state" and one can use a "small cage expan- 
sion" around typical liquid configurations in the spirit of 
the Einstein description of a crystal (in the DFT con- 
text, this would correspond to the self-consistent phonon 
approximation 3 ^). Truncated to quadratic (harmonic) 
order, this approximation leads to the following free en- 
ergy per particle 

3(1 — m) 

(p(m,A;T) = y — ' k B T [\n{2i:A) + 1] - mAC Uq (T/m) 

3 1 

- — k B T\n(m) H 4> Uq (T/m), 

Am m 

where A is the cage size, whereas Cu q (T) and 4>u q (T) are 
the expectation value of the Laplacian of the potential 
and the free-energy density for the equilibrium liquid (i.e. 
with m = 1) at temperature T. The free energy in the 
above equation should be minimized with respect to A (in 
the glass phase which we shall not consider here one also 
has to minimize with respect to m). The complexity (or 
configurational entropy per particle) in the equilibrium 
liquid is given by 
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(m,A*(m;T);T) 
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where A*(m;T) is the cage size at the minimum of 
Eq. ([3]). Because we only study the the temperature 
range above the ideal glass transition where m = 1, an 
expression of the complexity in this range is easily ob- 
tained by expanding Eq. ([5]) around m = 1, such that 
the definition in Eq. (|6]) yields 



£(T) = S liq (T) - - 



1 + ln 
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2(AV(r)>(T), 
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where (AV(r)) — 2np drAV(r)g(r) and A is the 
Laplacian. This expression shows that the complexity 
is obtained by substracting from the total entropy a "vi- 
brational" term requiring the sole knowledge of the pair 
correlation function g{r). 

Whereas the two approaches give in principle the same 
information, namely the complexity, the DFT is more de- 
manding as it requires a choice of trial configurations to 
describe the density profiles corresponding to the glassy 
metastable states. Except for some comments on the 
DFT, we shall therefore illustrate the ability of such cal- 
culations to reproduce the differing behavior of the WCA 
and LJ models by studying the complexity and Tk via 
the replica formalism and the small cage expansion. Wc 
only seek to give an idea of what the approach can or can- 
not do (with the commonly used approximations) and we 
therefore further simplify the computation by obtaining 
the g(r) for the WCA and LJ models via the HNC in- 
tegral equation for a one-component system. Of course, 
the choice of the HNC integral equation may have some 
quantitative impact on the results. 

We next consider an approach that mixes microscopic 
input and phenomenological considerations. Schweizer 
and coworkers 2 * have proposed a somewhat heuristic ex- 
tension of the MCT that incorporates activated mecha- 
nisms allowing ergodicity to be restored at temperatures 
below the T c of MCT. It focuses on single-particle dynam- 
ics and boils down to a stochastic nonlinear Langevin 
equation for the displacement of a particle r(t). The 
central quantity is the "nonequilibrium free-energy func- 
tional" which gives rise to the effective force exerted 
by the surrounding on a tagged particle. For a one- 
component liquid, this functional reads 



F eff (r)/(k B T) = 
d 3 q (S(q) 



-31n(r) 
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(2tt)3 p(S(q) + 1 



■ exp 



2 2 

q r 



(8) 



where S(q) is the static structure factor. In the viscous 
liquid regime, the competition between the two terms in 
Eq. ([5]) produces a minimum at small r and a maximum 
at a larger value. The relaxation time is given by the 
hopping time to escape from localization in the minimum 
and is expressed following the standard Kramers theory 



as 



t(T) = t exp 



Fb(T) 



knT 



(9) 



where Ff, is the height of the barrier between the maxi- 
mum and the minimum of F e ff(r) and the prefactor To 
includes, in addition to the short-time friction constant, 
information about the curvature around the minimum 
and the maximum. Since it has been shown to be weakly 
dependent on temperature and density, we take it as a 
constant. 

Finally, we also assess a correlation between dynam- 
ics and static pair density correlation that has recently 
been put forward, mostly on an empirical basis. Building 



on Rosenfeld's work 38 , Truskett and coworker a 22 i 39 have 
proposed a direct connection between thermodynamics 
and dynamics in liquids. In their picture, the transport 
and relaxation properties of a liquid are determined by 
the thermodynamic excess entropy (the entropy in ex- 
cess to that of the ideal gas, not to be confused with the 
"configurational entropy" in excess to that of the crys- 
tal). For simple atomic liquids in which only translational 
degrees of freedom are relevant, Truskett and coworkers 
went further to replace the excess entropy by its two-body 
contribution s 2 and therefore proposed a functional re- 
lationship between the diffusivity or the relaxation time 
and the static pair correlation function g(r) 22 ' 40 . For a 
one-component system, the two-body excess entropy is 
indeed defined as 



-s 2 (T)/fci 



cPr 



g(r)\n[g(r)] - g(r) + 1 



(10) 



which for a binary mixture is generalized to 
-s 2 (T)/k B = 



^2x a X/3 / d 3 r 



g a p{r) ln\g a /3(r)] - g a p{r) + 1 



where x a is the concentration of species a. We only focus 
here on the application of these ideas to the viscous liquid 
regime. 

The authors^ 2 - moreover suggested that a good empiri- 
cal expression to link transport and relaxation properties 
to two-body excess entropy is of the form 



t(T) oc exp[-Ks 2 (T)} 



(12) 



with K an adjustable parameter. (It should be stressed 
that the above equation is quite different than the much 
used Adam-Gibbs formula^, as the entropy comes in 
the numerator of the term in the exponential whereas 
the configurational entropy comes in the denominator in 
the Adam-Gibbs expression.) Note that this connection 
between excess entropy and relaxation has already been 
shown to fail in the case of associated liquids (square- well 
fluids^ 2 - and water models^ 3 -) and for confined liquids^ in 
the supercooled regime. 

All of the above descriptions of the dynamics of glass- 
forming liquids, at the level of approximation which is 
that commonly used but could in principle be improved 
(see below) , involve only the static pair correlation func- 
tions as nontrivial input, as shown by Eqs. ([21131), 
([11 [9]), and (110lll2p . In the following, we confront these 
approaches to the benchmark systems described above. 



IV. RESULTS 

A difficulty in assessing the ability of the microscopic 
approaches based on the knowledge of the static pair den- 
sity correlations to reproduce the differing behavior of 
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FIG. 1: Comparison of the relaxation times obtained by MCT 
(lines) and simulations (symbols) for the LJ (full lines) and 
the WCA (dashed lines) models at a density p = 1.2. 



the WCA and LJ glass-formers is that, when tested in 
simulation studies of glass-forming liquid models, these 
descriptions anyhow never provide a very accurate ac- 
count of the dynamical data. This is somewhat expected 
as the approximations involved in the theoretical formu- 
lations are quite uncontrolled and serious, while the phe- 
nomenology to be described is quite rich and complex. 
In approaches such as the MCT and the RFOT theory, 
a way to get around this is to implement some rescaling 
procedure to make predictions and simulation data on 
one type of observable fare as well as possible and use this 
rescaling to compare the predictions for other observables 
or phenomena. This is for instance what is usually done 
in test studies of the MCT: the critical temperature T c is 
rescaled in order to better match the simulation data on 
the relaxation time and comparisons are then made at 
relative distance from this temperature^. In the follow- 
ing, we shall allow for a possible rescaling of the theoret- 
ical predictions for one model, the LJ one, and use the 
same rescaling parameter for the WCA model. Doing 
otherwise by allowing for different rescaling parameters 
would completely side-step the issue of comparing these 
two systems, which have very close static structure fac- 
tors and very different dynamics at the same (p, T) state 
point. 

We first consider the results for the MCT. A detailed 
account has been given in Ref. and we illustrate here 
our findings for the sake of completeness. We have solved 
numerically Eqs. (J2J[3]) by using as input the partial static 
structure factors obtained in our simulation of the WCA 
and LJ models. From the solution for the intermediate 
scattering functions F a p (q, t) , we have extracted a relax- 
ation time r (defined from F aa (q,t — T)/S aa (q) = 1/e 
for the majority species a = A and for qoAA — 7.2) and 
we have determined the critical temperature T c . 

To begin with, we give in Fig. [T] an example of the dif- 
ficulty mentioned above: on an Arrhenius plot of the 
relaxation time at p = 1.2, the MCT predictions are 
quantitatively inaccurate even for the LJ model and one 



1 

• 


I 

rpL3 /rpWCA 


i 


i 

- 


■ 


/ J_jtJ / / V V V-' ji 












- 














-■--9 


-- ------e- * 


Replica / 








Mode-coupling | 


i 


i 


1 1.2 


1.4 

P 


1.6 


1.8 



FIG. 2: Comparison between the predictions of the MCT 
and of the replica-RFOT approach and the simulation results. 
Squares: ratio of the critical temperatures T c for the LJ and 
WCA models as predicted by MCT (open symbols) and es- 
timated from fits to the simulation data (closed symbols). 
Circles: ratio Tjf J /T^ CA predicted in the replica approach 
(open symbols) and of the ratio T ( ^ j /Tq VCA estimated from 
the Vogel-Fulcher-Tammann fit to the simulation data for the 
relaxation time (closed symbols). 

must rescale the temperature to obtain a more accept- 
able description. From the bare data, one can clearly see 
that the difference between the slowdown of relaxation of 
the WCA and the LJ models is considerably underesti- 
mated by the theory. To disentangle this from the overall 
quantitative inaccuracy, we plot in Fig. [2] the ratio of the 
critical temperatures T c of the LJ and the WCA models 
(which automatically accounts for a rescaling of tempera- 
ture that is the same for the two models) as predicted by 
the theory (see preceding section) and obtained from a fit 
to the simulation data in the range of temperature where 
this fit is possible^: the trends as function of density 
are similar but the theory is unable to capture the wide 
difference in the dynamics of the two models at typical 
liquid density (p < 1.4). 

For the replica (mean-field) approach of the RFOT the- 
ory, we have solved the HNC integral equation for the 
WCA and LJ models and used the resulting pair correla- 
tion function as input in Eqs. (J5l EI) to compute the com- 
plexity as a function of temperature and density as well 
as the associated Tk(p) at which the complexity goes to 
zero. The results showing the complexity as a function of 
temperature for several densities are given in Fig. [3] One 
observes that the curves for the WCA and LJ models are 
always very close (the difference is slightly bigger when 
one lowers the density to values close to the spinodal 
of the L J liquid) . To give some elements of comparison 
between theory and simulation, we have fitted the tem- 
perature dependence of the relaxation data of the WCA 
and LJ liquids with a Vogel-Fulcher-Tammann expres- 
sion, t = t cxp[BT /(T — T )]. As often done, we take 
To as an estimate for the Tk of the "entropy crisis" . We 
are not interested here in the absolute values but in the 
ratio of the temperatures for the WCA system and for the 
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LJ one. We plot in Fig. [2] the evolution with density of 
the theoretical prediction for the ratio T^ 3 /TW CA and of 
the simulation result for the ratio T LJ /T WCA . Without 
any doubt, the theory completely misses the widening 
gap between the two models as density decreases. 

Concerning the DFT approach to the RFOT theory, 
which we have not examined in detail, we just make the 
following remark. If, as done e.g. in Ref. 18(, the ap- 
proximation to the DFT approach relies on a WCA-like 
separation of the pair potentials with (i) the equilibrium 
structure (i.e., both the pair correlation function g(r) 
and the reference aperiodic lattice) determined by the 
truncated repulsive component and (ii) the attractive in- 
teraction treated in a mean-field-like fashion, then, es- 
sentially by construction, it cannot describe the differ- 
ence in behavior between the two models, with (LJ) and 
without (WCA) attractive tails. The fact that Hall and 
Wolynes 1 ^ use a somewhat different choice of cutoff than 
the conventional WCA one does not change this conclu- 
sion. The modified KRR truncation 4 ^ that they consider 
is equivalent to the WCA one for typical liquid densities 
(here p ~ 1.2); as implied from Fig. 5 of RefJ^, the 
predictions fare better at lower densities (corresponding 
here to p = 1.1) but they do worse at larger densities 
where they even fail to recover the merging with the LJ 
data found for the WCA truncation. 

To test the Schweizer-Saltzman approach for binary 
liquid mixtures, one should in principle consider a two- 
dimensional extension of Eq. (j8]) which takes into account 
the displacements of particles of both species. However, 
this is numerically very demanding. We have taken in- 
stead a shortcut that uses an effective one-component 
description focusing on the majority component A of the 
80:20 binary mixture. To do so we have implemented 
the procedure developed in Ref. 47 in a different context 
of a mixture of attractive and repulsive spheres. This 
approximation amounts to replacing p and S(q) in the 
one-component expression in Eq. ([5]) by pA and Saa(q), 
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FIG. 3: Temperature dependence of the predicted complexity 
(from the replica formalism and the small-cage expansion) of 
the LJ (full lines) and WCA (dashed lines) one-component 
models for several densities. A vanishing complexity defines 
the ideal glass transition temperatures, S(Tr-) = 0. 
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FIG. 4: Arrhenius plot of the relaxation time for the LJ and 
WCA models as predicted from the Schweizer-Saltzman ap- 
proach (open symbols) and obtained from simulation (filled 
symbols). For the theory, we show both the results of the 
effective one-component approximation with static structure 
data from simulation as input (lower curves) and those for 
the one-component models with pair correlation function ob- 
tained from the HNC integral equation (upper curves). Note 
that both theoretical predictions evolve with an unphysical 
sub- Arrhenius manner with temperature. 



respectively. Thus, we can directly feed the theory with 
the structure factors measured in the numerical simula- 
tions, much as we did to study mode-coupling theory. 
Additionally, we have solved the HNC integral equation 
for the one-component LJ and WCA potentials at the 
density p = 1.2, which we then directly use in Eq. ([8]). 

We plot in Fig.@]the predictions for the relaxation time 
as a function of 1/T at the density p — 1.2, together with 
the simulation data. The striking observation is that the 
theoretical predictions fail completely on two key aspects. 
First they fail to reproduce the super- Arrhenius behavior 
with its marked curvature. Second, they also fail to pre- 
dict any significant difference between the WCA and the 
LJ models. Puzzled by this major breakdown of the ap- 
proach, we have checked our implementation in the case 
of a hard-sphere system already studied by Schweizer and 
Saltzmar*2i. The results are presented in more detail in 
the Appendix and confirm the trends observed for the LJ 
system: the slowdown of relaxation as one cools or com- 
presses a glass-forming system is strongly underestimated 
by the Schweizer-Saltzman approach. It is true that the 
effective one-component treatment introduces an addi- 
tional approximation, but considering that we apply it 
to the majority component of the LJ and WCA mixtures 
which comprises 80% of the atoms, it seems more likely 
that the observed failure is intrinsic to the theoretical 
approach. 

Finally, we come to the correlation between the re- 
laxation time and the two-body part of the excess en- 
tropy. We have computed the latter from the simulation 
data for the partial pair correlation functions according 
to Eq. (jlip . The plot of r versus — S2//CB is shown in 
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FIG. 5: Relaxation time r versus two-body excess entropy 
— sijkB for the LJ (circles) and WCA (squares) liquids at 
density p — 1.2 (from simulation data). 

Fig. [5] for the typical liquid density p = 1.2. The data 
for the WCA and LJ models collapse at high tempera- 
ture (smallest values of — S2/^.b) but markedly deviate 
as T decreases (and —s^I^b increases), clearly showing 
the absence of master-curve. We note that the growing 
gap between the WCA and LJ data is somewhat reduced 
when plotted versus S2 in place of 1/T but the results 
nonetheless demonstrate that the two-body excess en- 
tropy does not uniquely determine the dynamics. Sim- 
ilar conclusions have been previously reached concern- 
ing supercooled associated liquid a 42 ' 43 and liquids under 
confinement 4 ^. 

We have also studied the evolution with density of the 
proposed correlation. It is expected, and indeed found, 
that the curves r versus — S2/&B for the WCA and LJ 
models converge as density becomes high enough. We 
have plotted in Fig.[6]the relaxation time (normalized by 
a high-temperature value r m which has a small residual 
dependence on density and temperature) versus —sij^-B 
for the two models. Besides showing the above men- 
tioned convergence at high density, one observes that the 
data for the LJ liquid at all densities collapse quite well 
on a master-curve, whereas a markedly different behav- 
ior is seen for the WCA liquid. This illustrates that the 
two-body excess entropy does not capture the difference 
in the density evolution of the slowdown of relaxation, 
which can be rescaled in the case of the LJ model but 
not for the WCA on o 23 ' 24 . This can be easily under- 
stood. Roughly speaking, and as far as fluctuations are 
concerned, the LJ liquid behaves as a soft-sphere model 48 
with power-law repulsive interactions characterized by 
an effective exponent 7. The relaxation time is then a 
function of the scaling variable p 1 /T only (as already 
sliow n 23 i 48 i 49 ), but this should also be true for the two- 
body excess entropy: one then expects a collapse of the 
r versus — S2/&B data, as indeed observed in Fig. [6] (this 
indicates that the adjustable parameter K in Eq. (TH21) is 
essentially independent of density). However, as the pair 
correlation functions of the WCA and L J models are very 



FIG. 6: Relaxation time r (normalized by a high-T value 
Too oc p 1 / 3 T~ 1/ ' 2 ) versus two-body excess entropy —s-zjkB for 
the LJ and WCA liquids at several densities (from simulation 
data). Note the good collapse of data for the LJ model (filled 
symbols) and its absence for the WCA one (open symbols). 

close, the two-body excess entropy of the WCA liquid has 
also an approximate scaling in p 1 /T. It is then obvious 
that the strong violation of density scaling which is ob- 
served in the slowdown of relaxation for this mode l 23 ! 24 
cannot be accounted for by the two-body excess entropy 
. This is what we see in Fig. [6] 

V. DISCUSSION 

The so-called "microscopic" approaches of the dynam- 
ics in glass-forming liquids that rely on the sole knowl- 
edge of the static pair density correlations fail when 
tested against the benchmark provided by two glass- 
forming liquid models with very similar static pair struc- 
ture and strongly different dynamics. This is so even if 
one adjustable rescaling parameter, the same for the two 
models, is introduced to improve the overall quantita- 
tive accuracy with respect to simulation data. Of course, 
by introducing more adjustable parameters, one would 
end up with a better description but this would obvi- 
ously be at the expense of the "microscopic" character 
of the approaches. Our demonstration is not flawless, 
as we have introduced in some cases an additional as- 
sumption in the form of an effective one-component de- 
scription of the binary mixture. However, we believe it is 
strong enough to seriously question the use of theoretical 
or empirical descriptions of the dynamics of glass- forming 
liquids that are based on the static pair density correla- 
tions only. This finding echoes with extreme theoretical 
models which have been devised where the slowdown of 
relaxation takes place in the complete absence of pair 
correlations among particles^— . 

Some of the approaches studied here could in principle 
incorporate more information on the structure than just 
the pair correlation function. This is not what is done in 
their common implementations, and this would certainly 
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increase the complexity of the computations while reduc- 
ing the simplicity of the message. Although we feel that 
merely introducing corrections due to the triplet correla- 
tion functions in the MCT or the DFT approach of the 
RFOT theory would not be sufficient, we cannot exclude 
that considering more involved approximations that ac- 
count for higher-order structural correlations would lead 
to a better description of our benchmark systems 5 ^. In 
the DFT approach of the RFOT theory, one could for 
instance envisage considering trial amorphous configura- 
tions that already account for the subtle structural dif- 
ferences between the WCA and the LJ liquid; or, in the 
sought correlation between dynamics and excess entropy, 
one could consider the exact excess entropy as in Ref. [39| 
in place of its two-body contribution. It remains to be 
seen whether this improves the predictions, and we hope 
that the present work will stimulate further investigations 
along these lines. 

We would like to stress again that the above conclu- 
sions apply to glass-forming liquids, for which it is well es- 
tablished that the dramatic slowdown of relaxation comes 
with only minor changes in the equilibrium pair density 
correlations. For fluids able to be taken at low temper- 
ature and low density, such as systems interacting with 
truncated repulsive potentials, possibly supplemented by 
short-ranged attractive interactions, a wide thermody- 
namic range can be covered within which the static pair 
correlation functions may show marked and nontrivial 
evolutio n 54 ' 55 . The ability of theories based on the static 
pair structure to grasp some qualitative and quantitative 
trends is then less challengeable in this case. 

Finally, we note that the present study has obviously 
nothing to say about theories of glass formation that rely, 
at least in principle, on more than the static pair density 
correlations, involving either higher-order static correla- 
tions or kinetic constraints (see Introduction). Such the- 
ories however are not microscopic in the sense that phe- 
nomenological input is necessary, either because the full 
blown theory is intractable (at present) or because some 
pieces connecting to the microscopic details are still miss- 
ing. A crisp test of the type proposed here is therefore 
harder to envisage. 



Appendix: Schweizer-Saltzman approach for hard 
spheres 

In this appendix, we check our implementation of 
the Schweizer-Saltzman approach for the binary LJ and 
WCA liquids by applying it to the hard-sphere model 
already extensively studied by these authors^. 

We obtain the static pair structure using three meth- 
ods: we solve both the HNC and PY integral equations 
for a monodisperse system of hard spheres to get S(q), 
or we follow the same procedure as in the main text and 
measure Saa{q) for the majority component of a binary 
mixture of hard spheres to use the same effective one- 
component approximation. We then inject this structural 
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FIG. 7: Top: Logarithmic plot of the relaxation time r versus 
packing fraction <f> for the hard-sphere systems. Predictions 
from the Schweizer-Saltzman approach with the g(r) obtained 
from the PY (squares) and the HNC (circles) closures for the 
monodisperse system, as well as from simulation of the binary 
mixture using the effective one-component treatment (trian- 
gles). We also display the simulation data (filled symbols) 
from Ref. [fi(J for the binary mixture. Bottom: same data 
plotted as a function of Z — P/(fceTp). 



information into Eq. © to obtain a prediction for the 
relaxation time which we represent either as a function 
of the volume fraction cf>, or as a function of the reduced 
pressure (or compressibility factor), Z = P/(ksTp). The 
latter is evaluated from the Carnahan-Stirling equation 
of state, which is quite accurate even in the high-(/> regime 
that we wish to discuss^. 

The results of the two theoretical approaches using 
HNC and PY are displayed in Fig. [7J and they fully 
match published results 2 - , which shows that we have cor- 
rectly implemented the Schweizer-Saltzman approach. It 
should be noted that when plotted in the representation 
log(r) versus Z, these data appear almost linear, which 
seems to suggest that in the regime relevant to most 
simulations and experiments, the Schweizer-Saltzman ap- 
proach seems to yield results similar to the free volume 
approach r ~ exp(Z), and that the asymptotic regime 
predicted by Schweizer— , r ~ exp(Z 2 ), is not entered. 

It is not possible to obtain numerical results for a 
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monodisperse system of hard spheres to assess these the- 
oretical predictions directly because the system crystal- 
lizes too easily. However, it is well-known that slightly 
polydisperse systems and binary mixtures can easily be 
compressed toward the glass transition with no sign of 
crystallisation and very little effect of polydispersity on 
the actual location of the glass "transition" 58 i 59 . Thus, 
we use numerical result s 60 ' 61 obtained for a 50:50 mixture 
of hard spheres with diameter ratio 1.4, and present them 
along with the theoretical results in Fig.[7J This compari- 
son appears reasonable as Schweizer and Saltzman them- 
selves compare their predictions with results obtained 
from simulations and experiments performed on polydis- 
perse systems. The comparison with numerically mea- 
sured relaxation times confirms the conclusions drawn 
from the study of the LJ potential that the Schweizer- 
Saltzman expression of the free-energy barrier consider- 
ably underestimates the slowing down of the dynamics. 
This is more easily seen in the bottom panel which shows 
that numerical data have a much larger "kinetic fragility" 
(they are well described 6 - 1 - by r ~ exp(Z a ) with a w 6) 
than the theoretical ones, as also found for the LJ/WCA 
systems. We note that the use of integral equations as in- 
put for the theory has two conflicting effects: on the one 
hand, integral equations tend to overestimate the struc- 
ture of the fluid (which should yield larger barriers), but 
on the other hand they are not sensitive to the proxim- 
ity of the jamming transition, a transition which is not 
captured by simple liquid-state approaches (close to jam- 
ming this would imply that barriers are underestimated) . 

Finally, as for the LJ and WCA models in the main 
text, we use the numerically measured structure factor 
Saa(q) for the main component (larger particles) and 
the effective one-component approximation for the free 
energy barrier. Although the approximation is expected 
to be less accurate in this 50:50 mixture than for the 
80:20 LJ/WCA systems, we find a trend very similar to 
the one reported in the main text, namely that integral 
equations seem to largely overestimate the evolution of 
the structure and that direct use of the "exact" struc- 
tural information yields theoretical predictions in even 
stronger disagreement with the simulations. 

This short study of the hard-sphere systems therefore 
validates our implementation of the Schweizer-Saltzman 
approach for the binary LJ and WCA liquids in the 
main text, and, as a way of consequence, the conclusions 
we have drawn from Fig. [4] 



Addendum to the appendix 

If the assumptions we made are correct, the last con- 
clusions of the above appendix suggest that the agree- 
ment reported by Schweizer and Saltzman between their 
theoretical results and hard-sphere data is largely due a 
cancellation of errors between a dynamical theory which 
underestimates barriers, and the use of integral equations 
which promotes them. As a further check of our assump- 
tions, we have studied the dynamics of a four-dimensional 
hard-sphere system. Indeed, in 4 dimensions, a monodis- 
perse system can be used with no crystallization issues 6 ?. 
In this case we can directly measure the 'exact' structure 
factor in the simulations and use the Schweizer-Saltzman 
theory to predict the dynamics, which can be compared 
to numerical measurements. In this comparison, the free- 
energy barriers are then estimated without approxima- 
tions (whether due to polydispersity or to the computa- 
tion of the structure factor) and no assumption has to be 
made to implement the theory. We use the data obtained 
by Charbonneau and coworkers in Ref. [(32| to obtain the 
results shown in Fig. [5] The results show the same fea- 
tures as in three dimensions, which therefore confirms the 
validity of our previous conclusions. 
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FIG. 8: Comparison of the dynamics of the four-dimensional 
monodisperse hard-sphere fluid obtained directly from sim- 
ulations and obtained by using the Schweizer-Saltzman the- 
oretical approach with the 'exact' (numerically determined) 
structure factor as an input. As in three dimensions, we find 
that using the exact structure produces free-energy barriers 
that are too small, and the observed glassy dynamics is poorly 
reproduced by the theory. 
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